Spin susceptibility of neutron matter at zero temperature 
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Abstract 



The Auxiliary Field Diffusion Monte Carlo method is applied to compute 
the spin susceptibility and the compressibility of neutron matter at zero tem- 
perature. Results are given for realistic interactions which include both a 
two-body potential of the Argonne type and the Urbana IX three-body po- 
tential. Simulations have been carried out for about 60 neutrons. We find 
an overall reduction of the spin susceptibilty by about a factor of 3 with re- 
spect to the Pauli susceptibility for a wide range of densities. Results for 
the compressibility of neutron matter are also presented and compared with 
other available estimates obtained for semirealistic nucleon-nucleon interac- 
tions and with more traditional many-body techniques, like Brueckner's or 
Correlated Basis Function theories. 

PACS numbers: 26.60,+c, 26.50.+X, 05.10.Ln 
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In this paper we show that the strong correlations induced by realistic nucleon-nucleon 
interactions reduce, by up to about a factor of 3, the spin susceptibility, Xi °f degenerate 
neutron matter. This reduction may have important implications for problems of astrophys- 
ical interest, like for instance, neutrino scattering rates in dense matter and, more generally, 
the study of supernovae and proto-neutron stars [|1],|2| . 

Starting from the pioneering work by Sawyer ||, several calculations of the neutrino 
mean free path in uniform nuclear matter have been performed 0Q-§[ which show that the 
effects due to strong interactions are relevant, particularly in the spin-density channel which 
couples with the axial vector current. A sizable reduction of x> or ) equivalently, a large 
value of the G Landau parameter leads to an appreciable suppression of the Gamow-Teller 
transitions. 

The neutrino momenta and the momentum transfers in many applications are small com- 
pared with the neutron Fermi momentum, and both energy transfers and the temperature 
are small compared with the neutron Fermi Energy. Therefore, the Landau parameters of 
neutron matter at zero temperature are the main quantities to compute in order to evaluate 
the mean free path of a neutrino in dense matter. 

In this respect, ab initio calculations of the Landau parameters or related quantities, such 
as the compressibility /C, the effective mass m* or the spin susceptibility x> f° r degenerate 
neutron matter are extremely important. Such calculations can now be performed because of 
recent advances in many-body methods, particularly those based on quantum simulations, 
and because of the much improved knowledge of the nucleon-nucleon interaction. 

Previous evaluations of the Landau parameters Fi and Gi (with I < 1) for neutron 
matter were based either on Skyrme-type potential models or on microscopic calcula- 
tions performed with semirealistic bare interactions [0-0. The qualitative behaviour of the 
compressibility ratio is 



K F ~ l + F 



(1) 



where K,f the compressibility of the non interacting Fermi gas, is similar in most of the 
various calculations. On the contrary, the spin susceptibility ratio, which is approximately 
related to G by 



Xf 1 + G, 



(2) 



o 



may differ up to a factor of 3 at equilibrium density of nuclear matter p = 0.16 fm -3 , and 
even more at higher densities. 

The Skyrme models generally predict much smaller values of Go with respect to micro- 
scopic calculations, and become unstable to spin oscillations, driving toward ferromagnetic 
ground states 0, since 1 + G becomes negative for densities in the range (2 — 4)p . At 
small momentum transfer, N-N correlations enhance the spin-response approximately by a 
factor {x/xf) 2 - Therefore, Skyrme models predict neutrino mean free paths smaller than 
microscopic models. 

On the other hand, the existing microscopic calculations have been performed with old 
semirealistic interactions, like the Reid or the Bethe- Johnston potentials, and no three- 
body force. In addition, the many-body methods used until now may be questionable for 
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the convergence of the underlying perturbation theory as well as for the treatment of spin- 
dependent correlations. 

We report the results of quantum simulations of neutron matter for the old Reid-t>6 
potential [|l(J and modern realistic interactions, based upon the Argonne vi$ two-body 



potential plus the Urbana three-body potential, UIX, (11]], denoted hereafter as AU18. It is 
well known that such interactions provide a realistic description of light nuclei and nuclear 
matter 0,0]. 

The quantum simulations have been carried out by using a new Diffusion Monte Carlo 
method, the so called AFDMC method [13 1. It introduces auxiliary field variables in the 



simulation, which makes it the only existing quantum Monte Carlo method that can handle 
spin-dependent nuclear Hamiltonians and a relatively large number of nucleons. In this 
approach the scalar parts of the Hamiltonian are propagated as in standard Diffusion Monte 
Carlo (DMC) [15|]. Auxiliary fields are introduced to replace the spin-isospin dependent 
interactions between pairs of particles with interactions between particles and auxiliary 
fields. Integrating over the auxiliary fields reproduces the original spin-isospin dependent 
interaction. The method consists of a Monte Carlo sampling of the auxiliary fields and then 
propagating the spin-isospin variables at the sampled values of the auxiliary fields. This 
propagation results in a rotation of each particle's spin-isospin spinor. 

The guiding function, \I/<r, in our AFDMC calculation is a simple trial function given 
by a Slater determinant of one-body space-spin orbitals multiplied by a central Jastrow 
correlation. The orbitals are plane waves that fit in the box times two component spinors, 
corresponding to neutron-up and neutron-down states. The overlap of a walker with this 
wave function is the determinant of the space-spin orbitals, evaluated at the walker position 
and spinor for each particle, and multiplied by the scalar Jastrow product. Such an overlap 
is complex, so the usual fermion sign problem becomes a phase problem. We constrain the 
path of the walkers to a region where the real part of the overlap with our trial function is 
positive. For spin-independent potentials this reduces to the fixed-node approximation. 

This method has already been applied to unpolarized neutron matter and neutron drops 
(A = 7, 8) with fairly realistic interactions that include tensor, spin-orbit and three-body 
terms. The neutron matter calculations have been done with up to 66 neutrons in a periodic 
box with a low variance (< 0.1 MeV per nucleon). The calculation scales in particle number 
roughly like fermion Monte Carlo with central forces |JT6|JT7[1 . 

We compute the spin susceptibility by applying a magnetic field to the system. Ignoring 
any orbital effects, the Hamiltonian is given by 

H = H -Y,Ji-b , (3) 

i 

where b = /i B and n = 6.030774 x 10~ 18 MeV/Gauss, and the susceptibility is defined as 

2 d*E (b) 



X = ~P\i 



db 2 



(4) 

b=0 



where p is the number density and Eo(b) is the ground energy in field b. 

Let us use the Pauli expansion of the energy per particle as a function of the spin 
polarization p = —dE (b)/db\ b=0 
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E(p) = E(0)-bp + ^p 2 E"(0) 



(5) 



where the derivatives are with respect to the polarization. Minimizing E(p) with respect 
to p one gets the following result for the spin susceptibility 



X 



H 2 p 



E"(0) 



(6) 



For a noninteracting Fermi gas the spin susceptibility is xf — p^mkf /(h 2 n 2 ). AFDMC 
allows us to get the energy eigenvalue, E (J Z , b), for the interacting system in a field b for a 
state of a given spin asymmetry J z — — iVj_. 

Assuming that the energy and polarization are known in terms of J z , E"(0) in Eq. (||) 
can be obtained as a straightforward application of the chain rule 



E"(0) 



dp 



d 2 E dE 



dJl dJ z 



dp 
dT z 



-i -i 



d 2 p 
dJ? 



(7) 



Since we are calculating the lowest energy state, the derivative of the energy with respect 
to J 7 vanishes. Therefore this result reduces to 



E"(0) 



dp 
dX 



2 d 2 E 
dJl 



(8) 



Let us consider the non interacting finite systems as a guide for the quantum simulations 
of the interacting ones. For such systems the energy is not a quadratic function of the 
external field b. In figure [j] we plot E (J z ,b), as a function of b, for four different systems 
with a finite number iV ~ 60 of noninteracting neutrons in a periodic box at p = 2pQ. The 
cases, for which we have done simulations, are shown in the figure, namely (iVj, iV|)=(33, 33), 
(33, 27), (57, 7) and (57, 0). One can see that the various E(J Z , b) are linear in b and each of 
them is tangent to the Pauli parabola (which refers to infinite Fermi gas case) at some value 
&o of the field b. We also compare in the figure the Pauli parabola with the exact result for 
the Fermi gas at that density. They are very close up to b ~ 50 MeV. Using Eq. (§), one 
gets x/Xf ~ 1 for J z = 50 and J z = 57. 

In the interacting case, the derivatives in Eq. (g) can be easily estimated by computing 
Eo(J z ,b) with AFDMC, and using the following equations 



dp E (J Z = J z0 , b = 0) - E (J Z = J z0 , b = b 



dJ x 



Jzo bo 



d 2 E _ E (J z = J z0 ,b = 0) - E {J Z = 0,b = 0) 



j 2 



(9) 
(10) 



r-0 



whose validity relies on the following reasonable assumptions: (i) for 6 = 0, Eo(J z ,b) is 
quadratic in J z ; (ii) for a fixed J z , E Q (J z ,b) is linear in b; (iii) the polarization is linear in 
J z . These assumptions become exact in the limit of an infinite system with J z and b small. 

The non-interacting case indicates the use of J z q = 50 and the value for 6 a ^ which 
E(J z q = 50, b) is tangent to the Pauli parabola (for the density p = 2p of Fig. [j], b = 53 
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MeV). Most of our calculations have been carried out with J z q = 50. At p = 1.25po we 
have verified the linearity of E(J z0 = 50, b) on b beyond b Q (1.25p ) = 39MeV finding that 
the result for x/Xf is largely independent on the value of bo. We have also verified the 
dependence of x/Xf on J z0 , performing simulations with J z0 = 6 and found that it is very 
weak. 

A time step At = 5 x lO^MeV" 1 was sufficient in most of the quantum simulations 
to obtain agreement between the mixed and the growth energies \TE\ within the statistical 
accuracy. 

We have made simulations with the Reid Vq interaction (Reid6), the same used in the 
Correlated Basis Function calculation of Ref. ||. Our result for x/Xf, at p = 1.25p , is about 
25% smaller, which indicates that the CBF perturbative calculations of Ref. |§, based on 
Jastrow-correlated basis functions, either have not reached a satisfactory convergence for the 
Landau parameter G or Eq (H) may not be sufficiently adequate for the Reid6 interaction. 

We have also considered realistic interactions, characterized by the so called v' 6 or v' s 
two-body potentials plus UIX three-body potential, and hereafter denoted as AU6' and 
AU8' respectively. The v' 8 is a two-body potential of the Argonne type which includes the 
four central spin-isospin components, plus the four tensor and spin-orbit ones. It fits the 
nucleon-nucleon experimental data and embodies the main features of the Argonne v is . 
The v'q potential is the v' s with its spin-orbit components removed. 

Our results for x/Xf-, obtained with the Reid6, AU6' and AU8' interactions are given 
in Table | and compared with previous microscopic calcuations of the same quantity. The 
Brueckner theory calculations of Ref. |J are performed with the full Reid potential, and, 
therefore are not directly comparable with our Reid6 calculations. However, we do not 
find a sizable contribution to x/Xf coming from the spin-orbit component of the two-body 
potential. Our results for AU6' and AU8' coincide within the statistical error. To address 
the problem of the influence of the three-body force on the results we have calculated the 
spin susceptibility at p = 2p using the Argonne v' & potential and no three body interactions. 
We have found that while the energy is reduced by roughly 25 MeV per particle, the spin 
susceptibility ratio is practically unchanged (0.31(1)). 

We have also calculated the compressibility /C, given by 

K 9 dp" + 9 dp ' 1 j 

where E (p) is a polynomial fit to the AFDMC energies E(J Z = 0, b = 0). For a Fermi gas 
the compressibility is JCp = 9n 2 m/ (kjh 2 ). The AFDMC results for /C//Cp, obtained with the 
AU6' interaction are shown in Table |J where they are also compared with the corresponding 
CBF estimates (AU6'-CBF) and other existing microscopic calculations fU^,[n|. 

For the sake of completeness we show in Fig. [| the equations of state of neutron matter, 
which have been used to compute the compressibility ratio of Table [El]. 

The AU18 results are taken from Ref. |13| , and have been obtained with the full AU18 



model interaction, by using variational FHNC/SOC methods. The AU6'-CBF [[19|] results 



have been obtained by using the AU6' interaction, as in the AFDMC simulations, and 



essentially the same many-body technique as in Ref. [13|. They also include the corrections 



coming from the lowest order elementary diagram, as discussed in Refs. [ PU| . |T7 |. 



5 



One can see that the AFDMC results for both the equation of state and the compress- 
ibility of the AU6' model of neutron matter are in reasonably good agreement with those 
obtained by using CBF theory. 

We have estimated the finite size effects in the AFDMC simulations by performing vari- 
ational calculations with the Periodic Box FHNC method of Ref. |[21|| . They indicate a 
correction which is at most 10% of the mixed energy per particle. An error of the same size 
is expected for the compressibilty ratio and of an order of magnitude smaller for the spin 
susceptibility. A more detailed discussion of the AFDMC calculation and the equation of 
state of neutron matter will be given elsewhere [|17[| . 

In conclusion, we have presented new results on the spin susceptibility and the com- 
pressibility of neutron matter at zero temperature, which show a strong reduction of these 
quantities with respect to their Fermi gas values. The calculations have been performed for 
realistic model interactions, which include tensor, spin-orbit and three-body terms, and by 
using the AFDMC method, a newly developed quantum simulation technique, capable of 
dealing with strongly spin-dependent interactions. Of particular relevance is the quenching 
by about a factor of 3 of the spin susceptibility with respect to the commonly used Pauli 
value, found at all the various densities considered. Such a reduction has strong effects on 
the mean free path of a neutrino in dense matter and should be seriously taken into account 
in the studies of supernovae and proto-neutron stars. 
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FIGURES 



FIG. 1. The energy of non interacting neutrons as a function of magnetic field b at p = 2po for 
various finite-sized closed shell trial functions with spin up and down values shown. Also plotted 
is the correct infinite system energy and the parabolic Pauli value. 
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FIG. 2. AFDMC equation of state of the AU6' model of neutron matter (dots); CBF the- 
ory [19] results for the same interaction model are in the shaded area where the highest values 
correspond to the variational estimate. The equation of state obtained in Ref.[13] for the AU18 
interaction by using FHNC/SOC theory is given by dashes. The errors are smaller than the sym- 
bols. 
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TABLES 



TABLE I. Spin susceptibility ratio x/xf of neutron matter. Our AFDMC results for the 
interactions AU6', AU8' and Reid6 are compared with those obtained from Refs. [8,9] by using 
Eq(|2|). The statistical error is given in parentheses. 



p/Po 


Reid § 


Reid6 § 


AU6' 


AU8' 


Reid6 


0.75 


0.45 


0.53 


0.40(1) 






1.25 


0.42 


0.50 


0.37(1) 


0.39(1) 


0.36(1) 


2.0 


0.39 


0.47 


0.33(1) 


0.35(1) 




2.5 


0.38 


0.44 


0.30(1) 







TABLE II. Compressibility ratio K./K-F of neutron matter. Our AFDMC results for the AU6' 
interaction are compared with those obtained with CBF theory [19] (AU6'-CBF) and those of Refs. 
[8,9,13] The statistical error is given in parentheses. 



P/Po 


Reid § 


Reid6 |] 


AU18 [ 


a 


AU6'-CBF 


AU6' 


0.75 


0.91 


2.06 


1.10 


0.85 


0.89(3) 


1.25 


0.70 


1.35 


0.71 


0.45 


0.47(3) 


2.0 


0.49 


0.77 


0.26 


0.23 


0.21(3) 


2.5 


0.42 


0.60 


0.15 


0.17 


0.14(3) 
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